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The invention relates to a method of analyzing a quantity having temporal and 
spatial variations, notably of analyzing perfusion of blood through the myocardivrai of a 
patient's heart 



Such a method of analyzing a quantity, notably of perfusion through the 
myocardium is known from the US-patent US 5 803 914. 



1 0 The known method displays information relating to the distribution of blood 

flow within the myocardium and information relating to the myocardixmi wall motion. This 
information on the myocardium is obtained in the known method from SPECT (single 
positron emission computed tomography) imaging data. 

In particular, by the known method a display ring of arc sections of the 

"15 myocardium is displayed with a color coding which depends on the perfusion ratio in each of 
the arc sections* The amount of the wall motion is represented by the thickness of the 
displayed arc. Inward and outward movement of the wall motion are mdicated by eTctending 
the vndih of die arc section at issue towards or away from the center of the ring resfpectively. 
The known method displays a series of several, e.g. five, images of sequential frames. Hence, 

20 it is cumbersome to gain a good rendition of the perfrision process through the myocardiiun. 
Notably, the user of the known method needs to examine each individual section of the 
myocardium in several images. Further, the number of sequential frames that can be shown is 
limited severely by the capacity of the monitor on which the frames are shown. 

An object of the invention is to provide a method of analyzing a varying 

25 quantity which furnishes an easy-reference rendition of the variations of the quantity. In 
particular, an object of the invention is to provide an easy-reference rendition of perfusion 
through the myocardium of a patient to be examined. 

This object is achieved by a method according to the invention wherein 
- a multidimensional output data array is formed 
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- lie multidimensional output data array comprises array positions arranged along at least a- 
first data-axis and a second data axis 

- values of the quantily are entered in the multidimensional output data array, such that 

- values of the quantity at substantially the same instant are entered at respective 

5 positions in the multidimensional output data array at equal positions along flie first 

data-axis and 

- values of the quantity at substantially the same spatial position are entered at 
respective positions in the multidimensional output data array at equal positions along 
the second data-axis, 

10 The values of the quantity at issue are dependent on at least two variables. 

According to the invention the values of the quantity at issue are arranged in the 
multidimensional ouQ)ut data array as a function of these respective variables along the first 
and second data axis respectively. For example, the quantity at issue may be dependent on 
spatial position and tune, a spatial position and temperature or the quantity at issue may 

is depend on a radial and a tangential position relative to a predetermined center. In particular, 
the invention is advantageously employed in that the values of the quantity at issue are 
arranged in the multidimensional output data array as a fimction of time along the first data- 
axis and as a function of spatial position along the second data axis. The multidimensional 
output data array can be displayed as an output image. When a two-dimensional ou^ut data 

. 20 array is employed, the two-dimensional output data array itself can be displayed as an output 
image which shows the content of the multidimensional output data array. When three or 
more dunensions are employed, for example two-dimensional cross-sections through the 
multidimensional output data array, or a rendering or two-dimensional projection of the 
multidimensional output data array can be displayed as an output image. In the output image 

25 tiie pixel-values or voxel-values like grey-values, color-values or brightness values, 

correspond to the values in the multidimensional output data array. The positions of the 
pixels or voxels in the output image correspond to the positions in the multidimensional 
output data array. 

The output image gives an easy-reference rendition of the time variation of the 
30 quantity at issue along the first data axis as the pixel variation in one direction in tiae output 
image. The output image also gives an easy-reference rendition of the spatial variation of the 
quantity at issue along the second data-axis as the pixel- variation in a second direction in the 
output image. Usually the first and second directions are orthogonal. Moreover, as the values 
of the quantity in neighboring spatial regions are displayed close to one another in the output 
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image, differences between the temporal variations of the quantity in respective spatial 
regions can be very conveniently monitored. In pardcular, according to the invention the 
easy-reference rendition requires very few images to be examined. When a two-dimensional 
output data array is employed, even a single output image is sufBcient to provide the easy- 
5 reference rendition of the variations of the quantity at issue. For example, when the invention 
is applied to perfusion through the myocardium, then differences in blood supply to 
respective sections of the myocardium served by individual coronary arteries are displayed 
in a single output image. 

These and other aspects of the invention will be further elaborated with 

1 0 reference to the embodiments defined in tiie dependent Claims. 

Preferably, for respective spatial sections, values for the quantity are acquired 
and entered into the multidimensional output data array at their appropriate positions along 
the second data axis. The values of the quantity for the respective spatial sections represent 
the spatial variation of the quantity. For example for individual spatial sections a local . 

15 average over the respective sections, a local minimum or local maximum value or the value 
at the center position of the spatial section at issue is used. It appears that in many 
applications such a representation on a per spatial section basis provides a fair representation 
of the spatial variations of the quantity at issue. Namely, the values of the quantity for 
individual sections give an accurate representation of the way the quantity varies spatially 

20 from section to section. The use of values for the quantity for respective spatial- sections 

achieves a reduction of the amount of data to be processed. Notably, when the values of the 
quantity are derived from a series of images, the amount of data involved in the 
multidimensional output data array is far less than the amount of data involved in the images 
themselves, especially when high-resolution images of a high di^nostic quality are 

25 * employed. The differences between the temporal variations of the quantity at individual 
spatial sections are easily viewed in the output image formed from the multidimensional 
output data array. Namely, the rendition of the values of the quantity along the first data axis 
forms time-traces of the quantity. Respective time-traces are formed for individual spatial 
sections and rendered at respective positions of the second data axis. 

30 . The amount ofdata is further reduced by using values ofthe quantity for 

respective time intervals. The value at respective successive time intervals provides a 
temporal sampling of the value ofthe for example continuously varying quantity. 

Preferably, values of the quantity for adjacent spatial sections are entered at 
adjacent positions in the multidimensional output data array and correspondingly shown at 
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adjacent positions in the rendition of the multidimensional output data array as an image. 
This allows spatially variations of the value of the quantity and of the temporal variation of 
the quantity to be clearly rendered visible in that these variations appear as relative sudden 
changes in the multidimensional output data array. 
5 In another preferred implementation values of the quantity for spatially 

radially adjacent positions are entered at adjacent array positions in the multidimensional 
output data array. Consequently, in the rendition of the multidimensional output data array as 
an image radial variations of the quantity are clearly displayed. For example when applied to 
the examination of perfusion through the patient's myocardium, radial variations of the 

10 perfusion are important indications for disfimctioning parts in the coronary arterial system. 

The values of the quantity at issue, notably perfusion through the patients 
myocardium are conveniently obtained from a series of diagnostic images. For example 
difiusion-weighted magnetic resonance images are used to obtain local values of perfusion of 
blood through tissue, such as the myocardium. From individual images, values of the quantity 

15 for several positions are obtained. From successive images, values of the quantity at 
successive instants are obtained, notably in that the successive images are acquired at 
successive instants in time. 

In a further preferred implementation of the method of the invention individual 
positions in the multidimensional output data array are linked to the corresponding positions 

20 in the images from which the values of the quantity are derived. Notably, a link is provided 
between the spatial section and the temporal instants from which respective values of the 
quantity are derived. Thus, when an interesting feature is noticed in the rendition of the 
midtidimensional output data array, the user may indicate that interesting feature. Then the 
method according to tihte invention automatically indicates tiie corresponding position in the 

25 corresponding image in the series. There are several ways to implement such a link. In a 
simple example the position of the indicated position on the displayed multidimensional 
output data array is calculated. From fliis calculated position in the multidimensional output 
data array the corresponding image in the series is derived as well as the corresponding 
spatial position in the image at issue, notably the position of the corresponding section in the 

30 image at issue. Then a marker is placed at tiie corresponding position in the image -at issue. In 
another implementation a look-up table is provided with represents the correspondence 
between positions in the multidimensional output data array and positions (of the sections) in 
the images in the series. The indication of positions in the displayed multidimensional output 
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data array can be done by way of a cursor in tiie form of an arrow or another symbol under 
control by a e.g. a moiise^ joy-sticky touch screen. 



multidimensional output data array the image at issue from the series is automatically 
5 displayed and the position is marked in the image linked to the indicated position in the 

multidimensional output data array. For a notable feature in the multidimensional output data 
array the corresponding position in the image at issue is then easily provided. On the other 
hand, a position in one of the images can be indicated and automatically the corresponding 
position in the miiltidimensional output data array is marked. For a suspect feature in one of 
10 the images at issue, the corresponding position in the multidimensional ou^ut data array is 
easily provided and the variation of the qxiantity at the suspect feature is easily examined. 



way to note intraesting perfusion features from the rendition of the multidimensional output 
data array and the method also provides an easy reference to the corresponding position in 

IS the image at issue that also shows the local anatomic details in which the noted perfusion 
feature occurs. This implementation is in particular advantageous when employed for • 
studying myocardial perfusion. Deviations from normal perfusion can be easily referenced to 
the relevant part of the anatomy of the myocardium. In particular, this implementation 
provides an advantageous tool for the cardiologist to detect local malfunctions in the 

20 coronary arterial system. . . , . 



system according to the invention is dejSned in Claim 1 1 . The data processing system is 
advantageously used to perform the method of analyzing the quantity such as perfusion 
through the patient's myocardium. For example the data processing system is coupled to or 

25 included in an magnetic resonance imaging system which acquires difiusion-weigfated 
magnetic resonance signals and derives the values of the quantity at issue. A computer 
program according to the invention is defined in Claim 12. The computer program according 
to tiie invention can be loaded into the working memory of a processor, such as a data 
processing system so as to enable the processor to perform the method according to flie 

30 invention. The computer program according to the invention may be made available from a 
data carrier such as a CD-rom or may be downloaded from a data network such as the world- 
wide web. 



Preferably, upon indication of a position in the rendition of the 



When applied to the examination of perfusion, the method provides an easy 



The invention also relates to a data processing system. A data processing 



These and other aspects of the invention will be elucidated with reference to 
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ibs embodiments described hereinafter and wifli reference to the aocompanymg drawing 
wherein 

Figure 1 shows a schematic representation of a cardiology image management 
system in which the method of the invention is employed. 
5 Figure 2 shows a diagrammatic representation of an exanq>le of the relation of 

a series of images to the multidimensional output data array formed according to the 
invention. 



I Q Figure 1 shows a schematic representation of a cardiology image management 

system in v\^ch the method of the invention is employed. The cardiology image management 
system includes an imaging modaUty 1 (IM) which suppUes image data to the data processing 
system 2 (DSP). The date processing system derives ihe multidimensional output data array 
from the image data. The multidimensional output data array is then applied to a display 

15 system. Also the image data from the imagmg modaUty 1 are supplied to the data processing 
system. The multidimensional output date array is displayed on the display system together 
with the image data. Optionally, image processmg may be applied to these image date by the 
date processing unit to improve the rendition of the image date on the display system. In the 
example shown in Figure 1, the imaging modality is a magnetic resonance imaging system 

20 which generates the image date in the form of perfusion images of for example the patient's 
myocardium. On the display system the series of myocardial perfusion images are shown 
with the multidimensional output date array. To this end the display system may comprise 
several display screens, or several display windows 31,32 may be opened on a display area. 

Figure 2 shows a diagrammatic representetion of an example of the relation of 

25 a series of images from to the muMdhneasional output date array formed according to the 
invention. 

In particular. Figure 2 shows by way of example one of the myocardial 
perftision unages 4. The image of the myocardium is divided mto several rings and individual 
rings are spUt into several sectors. The rings have the middle of the image of the left ventricle 
30 as then: common center. The sectors extend mainly tangentially, that is transverse to the 
radial direction from the common center. The sectors are numbered essentially spirally 
outward, that is in an mdividual ring, the sectors are numbered for example increasing in the 
clockwise direction and the numbering also increases when going from an inner ring to a 
more outward ring. Although only a single myocardial perftision image is shown in Figure 2, 
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ibis single image is one in a temporal series of myocardial perfusion images which show the 
temporal progress of the perfusion of blood through the patient's myocardium. Figure 2 also 
shows the multidimensional output data array S, corresponding to the series of myocardial 
perfusion images. In the example of Figure 2, the first data axis is the time axis (t) and the 
5 second data axis represents tiie number of the sector (#S). In the multidimensional output 
data array, the data values of the series of myocardial perfusion images are plotted. In 
particular, for respective sectors in the individual myocardial perfusion images the average 
brightness value is computed in the data processing system. These average brightness values 
present the local amount of perfusion of blood through the myocardium. In the 

10 multidimensional output data array as displayed, along the time axis, time traces of the 

perfusion process for the respective sectors are produced. Along the vertical direction, i-e, the 
axis representing the sector nimiber spatial differences of the perfusion in adjacent sectors are 
easily referenced. 

The positions in the multidimensional ou^ut data array are linked to 

IS respective sectors in the resi)ective image of the series. The user may indicate a specific 

position, e.g. where a striking change of bri^tness variations in the multidimensional output 
data array occurs. Such an indication may be done for example by a simple mouse-click at 
the position in the displayed multidimensional output data array. Upon indication by the user 
of a position in the displayed multidimensional ou^ut data array ,the myocardial perfusion 

20 image linked to the indicated position is put on the display screen and the sector at issue is 

indicated for example by the white arrow shown in Figure 2. The myocardial perfusion image 
linked to the indicated position may be displayed enlarged in stead of the complete series, or 
may be displayed in an additional wdndow opened on the display screen. The myocardial 
perfusion image linked to the indicated position provides additional information on the 

25 patient's anatomy, notably the structure of the tissue in the sector at issue in the area where 
the interesting perfusion event shows up in the multidimensional output data array. 
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1 . A method of analyzing a quantity having temporal and spatial variations, 
wherein 

- a multidimensional output data array is formed 

- the multidimensional ou^ut data array comprises array positions arranged along at least a 
S first data-axis and a second data-axis, 

- values of the quantity are entered in the multidimensional output data array, such that 
- values of the quantity at substantially the same instant are entered at respective 

positions in the multidimensional output data array at equal positions along the first 
data-axis and 

10 - values ofthe quantity at substantially the same spatial position are entered at 

respective positions in the midtidimensional output data array at equal positions along 
the second data-axis. 

2. A method as claimed in Claim 1 , wherein 

15 - values ofthe quantity are acquired for respective temporal instants and for respective 
spatial sections and 

- values ofthe quantity for individual spatial sections are entered at respective positions in 
the multidimensional ou^ut data array at equal positions along the second data-axis. 



20 3 . A method as claimed in Claim 1 or 2, wherein 

values ofthe quantity are acquired for respective time intervals and for respective spatial 
positions and values ofthe quantity for individual time interval are entered at respective 
positions in the multidimensional output data array at equal positions along the first data- 
axis. 

25 

4, A method as claimed in Claim I , wherein 

- values ofthe quantity for successive time intervals are entered at adjacent positions in the 
multidimensional output data array and 
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- values of the quantity for adjacent spatial sections are entered at adjacent positions in the 
multidimensional output data array. 

5. A method as claimed in Claim 4, wherein values pf the quantity for radially 

5 contiguous spatial sections are entered at contiguous positions in the multidimensional output 
data array. 

6. A method as claimed in Claim 1, wherein the values of the quantity are 
derived from a sesries of images. 

10 

7. A method as claimed in Claim 6, wherein values of the quantity at respective 
instants are derived from respective images in said series of images. 

8. A method as claimed in Claim 7, wherein respective positions in the 

1 5 multidimensional output data array are linked to respective spatial sections in respective 
images of the series. 

9. A method as claimed in Claim 8, wherein the miiltidimensional output data 
array is displayed, 

20 - a position m the displayed multidimensional output data array is indicated and 

- on the basis of the indicated position in the displayed multidimensional output data array 
the corresponding image of the series is displayed and the corresponding spatial section 
in the image is marked. 

25 1 0. A method as claimed in Claim 1, wherein the quantity pertains to perfiision of 

the myocardium. 

11. A data processing system arranged to- 

- form a multidimensional output data array 

30 - the multidimensioiial output data anray comprising array positions arranged 

a first data-axis and a second data axis 

- enter values of the quantity in the multidimensional output data array, such that 
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- values of the quantity at substantially the same instant are entered at respective 
positions in the multidimensional output data array at equal positions along the first 
data-axis and 

- values of the quantity at substantially the same spatial position are entered at 
respective positions in the multidimensional output data array at equal positions along 
the second data-axis. 



12. A computer program comprising instruction to 

- form a multidimensional output data array 



10 - the multidimensional output data array comprising array positions arranged along at least 
a first data-axis and a second data axis 

- enter values of the quantity in the multidimensional output data array, such fliat 

- values of flie quantity at substantially the same mstant are entered at respective positions 
in the multidimensional output data array at equal positions along the first data-axis and 

15 - values of the quantity at substantially tihe same spatial position are entered at respective 
positions in the multidimensional output data array at equal positions along the seppnd 



data-axis. 
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10 



In a method of analyang a quantity having temporal and spatial variations a 
multidimensional output data array is formed. THe multidimensional output data array 
comprises array positions arranged along at leastafirstdata-axisandaseconddata-ax^^ 

as a spatial axis and atemporal axis. Values of the quantity are entered in tiie 
multidimensional output data array. Values of the quantity at substantially tixe same instant 
are entered at respective positions in the multidimensional output data array at equal 
positionsalongtiiefirstdata^axis. ValuesofthequantityatsubstantiaUytiiesamespahal . 
position are entered at respective positions intixe multidimensional output data array at equal 

positions along the second data-axis. 



Fig 2 
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